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Abstract 



We extend the Bismut-Elworthy-Li formula to non-degenerate jump 
fT^ ' diffusions and "payoff" functions depending on the process at mul- 

^ I tiple future times. In the spirit of Fournie et al [14] and Davis and 

Johansson [10] this can improve Monte Carlo numerics for stochastic 

volatility models with jumps. To this end one needs so-called Malli- 

avin weights and we give explicit formulae valid in presence of jumps: 

(a) In a non-degenerate situation, the extended BEL formula represents 

possible Malliavin weights as Ito integrals with explicit integrands; (b) 

(*X , in a hypoelliptic setting we review work of Arnaudon and Thalmaier 

"^ ' [1] and also find explicit weights, now involving the Malliavin covari- 

1 -^ I ance matrix, but still straight-forward to implement. (This is in con- 

j^ . trast to recent work by Forster, Liitkebohmert and Teichmann where 

C^ ' weights are constructed as anticipating Skorohod integrals.) We give 

' some financial examples covered by (b) but note that most practical 

cases of poor Monte Carlo performance. Digital Cliquet contracts for 

instance, can be dealt with by the extended BEL formula and hence 

without any reliance on Malliavin calculus at all. We then discuss 

jS ' some of the approximations, often ignored in the literature, needed 

to justify the use of the Malliavin weights in the context of standard 
jump diffusion models. Finally, as all this is meant to improve numer- 
ics, we give some numerical results with focus on Cliquets under the 
Heston model with jumps. 



1 Introduction 

Modern arbitrage theory reduces the pricing of (non-American) options 
to the computation of an expectation under a risk neutral measure. It is 
common practice to assume that the risk neutral measure is induced by a 
parametric family of jump diffusions which can then be calibrated to liq- 
uid option prices. We can therefore assume that all expectations are with 
respect to a fixed pricing measure. A typical option on some underlying 
(St) then has (undiscounted) price 

E[f{ST„ST,,...,STj]=nfiS)]- 

For hedging and risk-management purposes it is crucial to understand the 
dependence on So and other model parameters. Computing 



a|E|/(5)i = E 



^^<^)f 



via finite differences can present computational challenges in Monte Carlo; 
just think of an at-the-money digital option near expiration. Broadie and 
Glasserman [7] showed that this problem is overcome by 

^IE[/(^)]=E[/(^)^] (1) 

where tt is the logarithmic derivative of the joint density of the random 
vector S_. On the other hand, the random weight it adds noise itself and it 
is important to localise: for instance by using (1) for an irregular, but com- 
pactly support and bounded, / and the usual finite difference technique 
for / — /, assumed to be nice (C^ will usually suffice). 

In two seminal papers, Fournie et al [14] and [15] use Malliavin calcu- 
lus to compute tt when no explicit transition density is known. They work 
with non-degenerate (or: elliptic) continuous diffusions but also cover some 
hypoelliptic situations. As is well known, elliptic results can be obtained 
by the Bismut-Elworthy-Li formula (Elworthy and Li [11], Bismut [6] ) and 
there are, in fact, other ways to obtain such results without Malliavin cal- 
culus: we mention in particular the idea of Thalmaier [27] of differentia- 
tion at the level of local martingales w^hich was employed by Gobet and 
Munos [18] in the present context. The point was that in many cases of 
practical interest, at least in absence of jumps, one does not need Malliavin 
calculus. (Specialists will note that Malliavin techniques are more flexible 
in the sense that different perturbations of Brownian motion yield different 
weights and there is an apriori interest to pick weights with small variance. 
In reality, it is hard to justify much effort in this direction as the potential 
gains are negligible to the improvements obtained by localisation.) 

Over the last decade it has become clear that pure diffusion models are 
unable to fit the short-dated smile and jumps have been included to mod- 
els to rectify this situation; Cont and Tankov [9] and Gatheral [16] provide 



two excellent accounts. The question has arisen as to how the above ideas 
can be adapted to models based on jump diffusion processes and we shall 
propose a quite simple solution to this along the ideas of Elworthy-Li by- 
passing both classical Malliavin techniques and its extensions to Levy pro- 
cesses that have been used in this financial Monte Carlo context. We note 
that a similar extension of the BEL formula, slightly less general than ours, 
was used recently by Priola and Zabczyk [22] to establish Lioville theorems 
for non-local operators. 

Let us briefly mention that in some cases a random weight tt can be 
constructed by conditioning arguments. Consider for instance the trivial 
example Xt = z + Bt + Nt, where i? is a standard Brownian motion and 
N a Poisson process. Conditional on Nt, any function of Xt is a (different) 
function of z + Bt, a pure diffusion with no jumps, and since the associated 
random weight n is universal (i.e. do not depend on the particular pay- 
off function) this also solves the problem for the jump diffusion X. This 
kind of reasoning leads immediately to the class of "separable" jump dif- 
fusions, considered in Davis and Johansson [10] via Malliavin calculus for 
simple Levy processes. We shall omit a detailed discussion since a refined, 
iterated conditioning argument can be used assuming only finite activity of 
the jumps (and without assuming separability in the sense of [10]). To this 
end, we quickly recall the BEL for continuous diffusions (see Section 3 for 
notation and assumptions) 



AE[/(xf,)]=E 



fix^r) I a{t)[R{t,xl)^] dWt 






where /^ a (t) dt = 1 and xq = z. Let < 5 < T be deterministic. The 
standard choice a = 1/T gives a weight, say ttq.t- Another weight (of 
higher variance) comes from a = 1/5 on [0, S*] , otherwise, and we call 
it TTo.s. One can also condition on xs and apply the BEL formula over the 
time interval [S, T], this yields another weight tts^t for the derivative of 
E[/(a:T)|a;5] w.r.t. xs- We leave it to the reader to check that, combined 
with the chain-rule, (d/dzj) = (dxs/dz) d/dxs, the weight ttq^t can be as- 
sembled from 7ro,s arid iis,T- In other words, instead of applying BEL on 
[0, T) one can apply it on [0, S) and [S, T). While we did not assume jumps 
in this discussion, it is clear that a cadlag discontinuity of x at time S does 
not pose a problem. This extends to any number of intervals and if we are 
dealing with a finite activity jump diffusion conditioning will reduce the 
problem to the one just discussed. The flaw with this sort of reasoning is that it 
makes fundamental use of a property which is completely irrelevant for the result 
to hold true: finite activity of jumps. In the general case, i.e. beyond finite ac- 
tivity, not only does the preceding argument break down, but jumps arise 
from a genuine stochastic integral w.r.t. a compensated Poisson random 
measure and any conditioning on jumps must fail. 



On the other hand, by maintaining a finite activity assumption on the 
jumps and some conditions on linkage operators, the eUipticity condition 
has been relaxed to hypoellipticity by Forster, Liitkebohmert and Teich- 
mann [13]. Unfortunately, their 'linkage' condition on the jump vector 
fields excludes many examples of financial interest B The main contri- 
bution of [13], in our view, is to establish new conditions for integrability 
of the inverse of the Malliavin covariance matrix C in presence of jumps. 
Recent progress in this direction was also made by Takeuchi [26] who man- 
ages to bypass Norris' lemma, which, in a sense is the bottleneck of the 
arguments in [13]. Thus, noting that criteria for integrability of C^^ are 
available in the literature, and can also be checked by hand in many exam- 
ples, we show that suitable integrability of C^^ allows to extend a recent 
result by Arnaudon and Thalmaier [1] and we so obtain non-anticipating 
Malliavin weights for hypoelliptic diffusions with jumps of possibly infi- 
nite activity also allowing for the 'linkage' condition in [13] to be relaxed. 

It is worthwhile to ponder for a moment which financial examples re- 
ally benefit from BEL / Fourni et al type formulae. The standard hypoel- 
liptic example in finance is an Asian option but computation of Greeks 
with (intelligently chosen) finite difference perform rather well. In fact, 
most jump diffusion models used in practice have an essentialljo elliptic 
diffusion part and also quasi-closed form expressions for European option 
prices and the usual Greeks, typically by Fourier methods i.e. by low- 
dimensional integration. Thus, the focus should really be on instruments 
without (quasi-)closed form prices for which finite difference methods per- 
form poorly. In fact, there is a very popular family of such contracts in eq- 
uity markets, namely digital cliquets, and the numerical difficulties for risk 
management are well-known to practitioners. Surprisingly perhaps, there 
seems to be no result in the literature that applies to computing sensitiv- 
ities of digital cliquets under, for instance, the Heston-model with jumps: 
the separability conditions of [10] are far too stringent, the relevant state- 
ment in [13], Proposition 1 to be precise, still contains the (here unneces- 
sary) linkage condition which is not satisfied|jnor do Heston-type models 
satisfy the strong C°° assumptions of [13] Q 

This paper is organised as follows. We prove that the Bismut-Elworthy- 
Li formula holds for a generic non-degenerate time-inhomogeneous Marko- 
vian jump-diffusion; n is given explicitly as a stochastic integral involv- 



^Indeed it is easy to see tliat ttiis condition fails in the case where the jumps in the stock 
are log normal as in the Merton model (or any example in which the Levy measure has full 
support). 

^The gap between elliptic and what one has in some real examples is subject of Section 5 
of this paper. 

^One could re-run the MaUiavin calculus arguments of [13] in the elliptic setting to get rid 
of this condition or, in fact, make rigorous the iterated conditioning argument outlined above. 

^This is just to say, that approximation arguments similar to those discussed in Section 5 
of this paper would be needed. 



ing the flow and the right-inverse of the diffusion matrix, just as in the 
classical Bismut-Elworthy-Li formula (which is recovered in the absence 
of jumps). A similar presentation is given for the second derivative. In 
section 4 we demonstrate how Malliavin calculus may be used with an 
appropriate choice of perturbation to provide explicit weights in a hypoel- 
liptic setting. In section 5, as a case study, we show how to represent the 
spot sensitivit}0 in the Heston model with jumps (also known as SVJ) and 
the Matytsin double jump model (an extension of SVJ and also known as 
SVJJ). Both models are described in detail in Gatheral [16] and are popular 
in the industry because of their quasi-closed form solutions for European 
options in terms of Fourier-transforms, which we use for numerical bench- 
marks for some simulations in the last chapter. An honest application of 
the BEL formul4j to these (and many other practically relevant) examples 
requires approximation argument which, in our view, have been neglected 
in the literature. 

The authors would like to thank James Norris, Chris Rogers and Anton 
Thalmaier for related discussions. 



2 Preliminaries 

We collect some background material from Gikhman and Skorohod [17]. 
Our focus is on the strong solution of 



z+ I Z{s,xl_)ds+ / X{s,xl_)dWs 

. ^° (2) 

/ Y{t,xl_,y){n-iy){dy,ds) 
Je 

where Wt = (W^, . . . , W/") is an M^-valued Brownian motion on some 
probability space (r2,JFt,P) and /.t is a (r2,JFt, P)-Poisson random measure 
on _E X [0, oo) for some topological space E such that u, the compensator 
of fi, is of the form G{dy)dt for some cr-finite measure G. The vector fields 
Z{t,x) and Y{t,x,y) e R'^, X{t,x) e /:(R™,R'^) for all t G [0,r], x e 
W^ and y G E. We will always assume at least the following conditions 
which guarantee the existence and uniqueness of a solution to the SDE 
(see Gikhman and Skorohod [17]) 

1. For all X eM.'^ and t e [0, T] 

\Z{t,x)\^ + \X{t,x)\''+ [ \Y{t,x,z)\''Gidz)<Cil + \x\^) 
Je 



^Delta would be a misnomer here. Cent and Tankov [9] contains a nice discussion of how 
these concepts are related in presence of jumps. 
*... and Malliavin techniques in general ... 



2. For all x.zeR'^ and t G [0, T] 

\Z{t,x) - Z{t,z)f + \X{t,x)-X{t,z)f 

+ f \Yit, X, y) - Yit, z, y)l'G(dz) < C\x - z\^ 



Throughout we fix the option expiry time T > and consider a payoff 
/ : M** ^ K, and we frequently work with the process Xj '^ for t' < t 
defined as the solution to the SDE 



c*''^ = z+ / Z{s,x*^l')ds+ / X{s,x*^l')dWs 



f Jt 

t 



(3) 



Y{s,xll%y){fi~ iy){dy,ds). 



JE 



We will write / £ Cl;{W) to mean that the function f : W ^ R is k- 
times continuously differentiable with / and all its derivatives up to order 
k uniformly bounded. Our method of proof will rely on ensuring that the 
function u{t, z) : [0, T] x K'' -^ R given by u{t, z) = E[/(a;j.^)] satisfies the 
backward equation of Kolmogorov. 

dii 

gtu{t,z) + —{t,z)^o. 

when / e C^i'R'^). It is well-known that for / e Cl'^i[0, T] x R'^), f lies in 
the domain of G and 



+ nf{t,x + Yit,x,y))- fit,x)^Y.^\t,x,y)^^{t,x)\G{dy) 

(4) 

where A = X{t, x)X'^{t, x) (here X'^ denotes the transpose of the matrix 
X). The arguments leading to the proof of the following theorem may be 
found in Gikman and Skorohod [17]. 

Theorem 1. (Kolmogorov's backward equation) Let x\ '^ for t' < t < T repre- 
sent the solution to the SDE (3) with the vector fields Z, X and Y satisfying the 
existence and uniqueness conditions. Suppose further that for every t <E [0, T] and 
y & E the following conditions hold 

Z{t, .), Xit, .) and Y{t, .,y) e CU^") (5) 



The functions Z{., .), VZ(., .), y^Z{., .), X{., .),VX{., .), V^X{. 



and / \Y{.,.,y)\'G{dy), \^Y{., .,y)\'G{dy)J \^'Y{., .,y)\'G{dy) 
Je Je Je 

are all continuous on [0, T] x R'^ (6) 

sup f / (|y(i,x,2/)|'= + |Vy(i,x,2/)|'=)G(dy)) <oo (7) 

(t,a;)e[0,T]xR'' \J_B / 

/or k = 2, 3, 4. Then iff e C^{W^) the function u{t, z) = E[/(x^^)] is such that 
u G Cfa '^([0, T] X M'') and sah's/ies 

— {t,z) + gtu = Q 
with boundary condition linif^T u{t, z) ~ f{z) and where Qt is given by (4) . 

3 The Bismut-El worthy-Li formula for 
j ump-dif fusions 

The argument of Elworthy and Li [11] extends in a straight-forward way to 
jump-diffusions. 

Theorem 2. Fix some T > and consider x^ = x^'' the solution to SDE (2) 
on the interval [0, T] and suppose that the conditions of Theorem 1 are satisfied. 
Further assume that the diffusion matrix X{t, x) has a right inverse R{t, x), and 
satisfies the following uniform ellipticity condition 

y'^Xit,x)X^{t,x)y>e\y\^ 

for every t € [0, T], x,y € R'^ and some e > 0. Then, if a £ L^[0,T] is any 
deterministic function which satisfies 

T 

a{t)dt = 1 



'0 



and f e C^(R'') the following is true for alll < k < d 

riT.'r \ 

(8) 



Ae[/(x^)] = E 

OZk 



fix^j.) 1^ a{t)i^R{t,xt_)^j dWt 



Moreover, if we consider < Ti < . . . < r„ < T and a function of the form 
/(xf,^ , . . . , x'f^ ), where f e C^{R'^ x . . . x R"^) and let a € L^ [0, T] be a deter- 
ministic function satisfying 



Ti 

a{t)dt = 1. 





Then, for alll < k < d, the following is true 
d 



dzk 



E[/(xf,^,...,a;f,J] 



E 






dWt 



(9) 



Remark 1. In the absence of jumps and with a{t) —T^on [0, T] we recover the 
classical Bismut-Elworthy-Li formula. 

Remark 2. It is easy to see that the uniform ellipticity condition gives rise to the 
fact that Jg a{s) (r{s,x^_) q^^ ) dWs is a martingale on [0, T]. To see this 



take t = s, X = Xg^ and y = R{s, Xs- 






and observe that 



R{s,xU- 



dxl 



dzk 



< e 



dxl 



dzk 



a.s. 



Consequently, 



„ owl^ft-.'-)^ 



dt 



<e-^E 



sup 

0<t<T 



dx^ 



dzk 



a{t)'^dt < oo. 



Proof. For t < T we apply Ito's formula to the function 



for t < T. Since 



uit,z)^E[fix'/-)]^PT-tf{z) 



du , ^ 



the ds-term vanishes, leaving only a constant and the two martingale terms. 
Letting t ^ T we find 



/(x^) =u(0, x) + / (Vu(s, xUX{s, xUY dWs 
Jo 

/ (w(s,a;^_ +Y{s,xl_,y))~u{s,xl_)){n~v)ids,dy) 
'0 Je 

(10) 



The integral featuring above with respect to /i — i^ is a discontinuous L^- 
martingale which is orthogonal to the martingale J ais) ( R{s, xl_ ) -gj— ) dWg 



Multiplying by J^ a{s) (^R{s, xl_ ) ^ 
the result 



dWs and using Ito's isometry gives 



E 



/(•^t) 



a{s)(R{s,xl^)^'' 



dzk 



dW, 



Ox 



a(s)E 



Vu(s,xf_) 



dzk 
dxl_ 



_d_ 
dzi 



«[/(= 



■Ti7- 



■.^tJ] 



_9_ 

dzi 



ngi^k)] 



E 



.9(a^Ti; 



Ti 



dxf 



, «W^^(t,.?_)^ 



dW^t 



which concludes the proof. 



dzk 

^ d 

ais)— — Efw(s, xl_)]ds 
dzk 

a{s)^E[f{x'T)]ds 

Cl^k 

i-m^T)] (11) 



ds 



We justify the progression from the second to third line by a routine ar- 
gument based on the boundedness of Vu and V^u and the definition of 
-^ as the L^-limit (as h —> 0) oi the random variables h^^{xf'^ ^'^'' — Xf) 
for fixed t. Also, we justify the third to fourth line in (11) by the obser- 
vation that u{t,xl) -^ u{s,Xs^) almost surely as t f s and so bounded 
convergence gives E[u(t,a;j)] -^ E[u{s,Xs-)]- But for each t e [0,r]we 
haveE[u(t,a;^)] = E[E[/(xf,)|.Ft]] = E[/(a;^)], so E[u(s,x,,_)] = E[/(a;f,)]. 
For the final part, we note that the function 5 : M^ ^ K defined by g{x) = 
E[f{x, Xrj}^'^, . . . Xj}'^)] has the property that g e C^(E^) and, moreover, by 
the Markov property 

.g(xf,J = nii'^T,,- ■ . ,x?.J|a(xf,J] = E[/(xf.^, . . .,x'rJ\J^TA a.s. 
Consequently, by (8) we have 



D 



Remark 3. Under stronger condition on the vector fields ( see Theorem (2-28) 
of Bichteler, Jacod and Gravereaux [5] ) we can ensure the existence of a density 
Pt{z, y)for the random variable x^ with pt S C^{R'^ x M''). It is then possible 
to relax the regularity restrictions on f so that we need only make a measurability 
assumption on f. 

Remark 4. The result (9) extends the result of Section 3.2 Fournie et al [14] to 
jump-diffusions. We notice that the form of the results do not correspond exactly, 



their weight is represented by 

with a e L^[0,T] satisfying J^ ' a{t) ~ Ifor all 1 < i < n, and our weight n is 
a particidar case of this when a = on [Ti,T]. However, it is clear that ifa^O 
on [Ti,T] then 

Var^n) < Var^n). 

Since the efficiency of Monte Carlo is optimised by the choice of the minimal vari- 
ance weight we would always choose a = on [Ti , T] and hence there is no 
conceivable practical advantage to representing the weight by tt. 

We may adapt this approach to deal with higher order derivatives as 
well. 

Theorem 3. Suppose that XX^ is uniformly elliptic and further assume that the 
conditions on the vector fields are strengthened so that the following conditions 
are satisfied. For every t G [0, T] and y & E 



Z{t, .), Xit, .) andY{t, .,y) e C^iR""). 



(12) 



For every I eNU {0},S/^Z{., .),\/'X{., .) and [ \\/^Y{., .,y)\^G{dy) 

Je 

are continuous on[0, T] x W^ (13) 



ForreN with r >2 and I ^ 1,2. 

sup ( [ {\Y{t,x,yW + \y'Y{t,x,y)nG{dy))<^ (14) 

(t,a;)e[0,T]xR<' \J E / 

Then, if f e Cf (R'*) and, for each t G [0, T], R{t, .) G C^ (R"*) (where the bounds 
on R{t, .) and Vi?(i, .) hold uniformly in t g [0, T], the following formula holds 
for all 1 < j,k < d 



&" 



dzjdzk 
4 



E[/(^t)] 



y2 



E 



dxi 



T/2 \ ^-^j 



T/2 



dxf 



dw^l^ [Rit,.U-g^^ 



dWt 



-E 



^E 
T 



/-T/2 / 



dxf_ dx^_ 
dzj dzk 



dWt 



■T/2 / Q2 2 



dWt 



10 



Moreover, if we consider < Ti < . . . < r„ < T and a function of the form 
f{x^ , . . . ,x^ ), where / e C^(R'' x . . . x R''). Then, the above result remains 
true when we replace f{x^) by f{x^ , . . . ,x^ ) and T by Ti in the above formula. 

Remark 5. The conditions on the vector fields are stronger than needed, but we 
state them in their current form for simplicity. 

Proof. Define the function w : M'^ x M"^ — > M by w{x,y) = V f{x)y. Then 
it is easy to verify that the function p : [0, T] x R'* x K'' -> R defined by 

p{t, z,y) = E It; I x^^, gj j satisfies the backward equation associated 
to the generator of the R'' x R"^' -valued diffusion ( xl , -^ ) . The argu- 
ment now proceeds as before; applying Ito's formula to p ( i, Xj , -^ J , let- 

ting t ^ T and then multiplying by /„ I R{t, xl_)-^ J dWt and taking 

expectations allows the argument to be concluded as in Theorem 2.3 of 
ElworthyandLi [11]. D 



4 Relaxing the EUipticity Criterion 

For simplicity we now assume that the vector fields are time homogeneous. 
We denote by U the £(R'',R'')-valued process given by Ut = V^xf and 
denote its inverse, when it exists, by Vt- We define the Malliavin co variance 
matrix 

Ct{z) = / {Vs^X{xl_)){Vs-X{xl_)fds 







and make the following a standing assumption. 

Assumption 1. For fixed T > 0,Ct is invertible a.s. and moreover |C^^| e L^ 
for all p > 1. 

This assumption is known to be true in certain cases, for instance it 
holds in the diffusion case under Hormander conditions on the vector 
fields (see Nualart [21]), and more recently it has been shown to hold in the 
jump diffusion case for finite intensity jumps under uniform Hormander 
condition (see Forster, Liitkebohmert and Teichmann [13]). For more gen- 
eral jump processes the problem is more involved but ideas in this setting 
have been developed in Cass [8] and Takeuchi [26]. 

We now prove an extension of Theorem 3.2 of Arnaudon and Thal- 
maier [1] which allows us to give an explicit representation of the weight 
in terms of an adapted R'^-valued process . Note that the result of Forster, 
Liitkebohmert and Teichmann [13] where the weight is given in the form of 
a anticipating Skorokhod integral may be converted into sum of integrals 
of adapted processes using the expansion formula (1.49) in Nualart [21]. 
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A representation of this type is more desirable from the point of view of 
simulation. First we recall some concepts from Malliavin calculus. Let a be 



an 



-valued previsible process such that for T > fixed 



E 



exp I t: 



\ash\ ds 



< oo, ft, G M'' locally at 



(15) 



and define a new probability measure by 



Z'^ 



J=T 



cxp 



< ai,h,dWs > ~7: 
2 Jo 



\ash\ ds 



and let Zt ~ E[Z!^\J^t\ for < i < T. Introduce a perturbed Brownian mo- 
tion dW^'' = dWt + ashdt and denote'byXf,Cf{z) the processes correspond- 
ing to Xf and Ct{z) when the underlying Brownian motion is replaced by 
Wl^. The crucial ingredient to the following result is the observation that 
the change of measure decribed above has no effect on the Poisson random 
measure. 

Theorem 4. Suppose Assumption 1 is in force along with the following condi- 
tions on the vector fields 



and swp / \V"Y{x,y)\'^G{dy) < oo for all n eN. 

xGS.'' Je 



i'^) 



Further assume 



snp snp\{I + '\/xY{x,y)) ^|<oo. 

xeMiiyeE 



(16) 



Then, for any f <E C^ (M^) and j e {1, 2, . . . , d} zve have 



dz. 



■E[/(xf.)] = E 



I{^t) I 



Vt-X 



{xt-)dWtY C:^\z)ej 



-i:(-^-(=)(^ 



ci^{z))c^\z) 

h—O / 



k.jy 



Remark 6. Condition (16) is there to ensure both the existence ofV = U^^ and 
that V e LP for all p > 1. In practice this can often be relaxed in favour of some 
less stringent condition (see Example 2 below). 

Proof. The fact that /i remains a Poisson random measure with compen- 
sator 1/ under P'' follows from Theorems (3.15) and (3.34) of Jacod [20] . 
We then observe,since x is a strong solution to the SDE (2), that 



T — 



h=0 



E[/(x^)Z^(C^(z)-i)fc„ 



(17) 
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Choosing the perturbation 



with an increasing sequence of previsible stopping times (r,i) chosen such 
that a" satisfies condition (15) and such that r„ f T.An elementary appli- 
cation of Ito's formula can be used to show 

d 
dhk 



x^ = Ut Vs-X{xs-)iVs-Xix,.)y dsek = UrCr^ek- 

h=0 Jo 



Using this we may expand (16) to get 



E[/(2;f.)C/Ta„C^^ej]=E 



fi^r) ( ( / " Vt-Xixt-)dWtYc^\z)ej 



fc=i / 



We let n ^ oo and expand the second term on the right hand side to give 
the stated result. D 

Example 1. (Bachelier with jumps, Asian options) We assume 

dSt = (jdWt + dNt 
dAt = St-dt, 

with some Poisson process N of finite rate. The Malliavin covariance matrix has 
the particularly simple form 

and so the second term on the right hand side of the formula in Theorem 4 drops 
out leaving us with 



^niiSr.Ar)]^^^ 



j; I Wtdt 



\wt\ 






USt.At) 



Example 2. (Exponential Levy, Asian options) Consider the following model for 
the evolution of a stock price 



dSt = P St-dt + aSt-dWt + / ySt-iji - v){dy, dt) 

Jy>~l 

dAt = St-dt, 
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ivith Ao = 0, and where W is a Brownian motion and /i a Poisson random mea- 
sure with compensator ^{dy, dt) = G{dy)dt. Make the assumptions that for all 
p> I and arbitrary S > 

{l + y)~PG{dy)<oo and f {I + yfG{dy) < oo. (18) 

-1 Jy>l 

Condition (16) is not satisified in this example, however it is easy to show by 
truncating the jumps at some arbitrary level that the theorem may be applied. 
Assumption (18) may then be invoked to guarantee the resulting formula remains 
valid in the limit as the truncation parameter goes to zero. We notice also in this 
case that the vector fields are not bounded and similar approximation results are 
needed, details on how this type of argument can be made rigorous are given in 
the next section but we omit them here for the purpose of clear exposition. The 
Malliavin covariance matrix may be computed 



Ct — 



2q2 



-'T^SoAt 



-a^SoJoAtdt^ 






It is easy to show that 



and then 



d 

dhTi 



_d_ 

dhi 



h=0 



Sf = a SoStt 



dh-2 



Si' 



'O^St 



h=0 



A.,ds 



h=0 



Al 



(T^So J^ sSs-ds ifi ~ 1 

a^dl^Aids-Ajl^Ads) if I ^2 



We nohce that detCt ^ a^S^ I t /„ A^dt - ( /„ Atdt) j , and 

C^' ^ {detCr)-' ( f C^* ^' "'^"f, '^''^'] . 

We must show that Cf,^ <E L^ for allp> 1. To see this it suffices to check that 

P {detCr < e) is o(eP) as e -> 0. 
To this end we note that for any < S < 1 



' {detCr < e) < P detCr < e, inf St > S, sup St < S'^ 

0<t<T 0<t<T 



inf St<S 

0<t<T 



sup St > S 

0<t<T 
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We now show that on the set B :— {info<t<T St > 5, supo<t<7^ St < S^^} we 
have, for 5 ~ (5(e) appropriately chosen , detCx > £• To see this note that on B we 
have St < At{Lu) < S^^tfor all t e [0, T], and so we define 

A^{f : [0,T] ^M, such that 5t< f{t) < 5~^tforallte [0,T]}. 

Then by examining the form of the determinant we have 

detCr > a^S^T'^ inf varf{U) 

where U ^ Uniform]}^, T]. We may bound the left hand side from below by Cheby- 
shev's inequality, so that for any a > 

varf{U)>a'n\f{U)-nf{U)]\>a) 
>a'P(/(C/)<-a + IE[/(C/)]) 



and taking a = ^ gives varf{U) > ^P {f{U) < ^). Since f{t) < 5-H 
[0, T] we have P (/([/) < ^) > P (c/ < ^) which gives 

g4rp3 



in 



varf{U) > 



64 



andsodetCr > ^ '^J '^" :== CJ". Choosing 5 ^ {C-^e)^/'^ we see that detCr > 
e on B. It therefore suffices to show that 



inf 5** < e and^i sup St > e^' are o(eP) as e ^ 
o<t<T J \0<t<T J 

for every p > 1. We show this for the infimum, the supremum being a simple 
modification of this argument. To this end we write St — S^e^* where Xt is a 
Levy process with triplet {a'^jfl, G) with 

il = ^L~]-a'^~ [ yG{dy), G{A) = G^e^ ~ 1 : x £ A}) for A e BiA). 
2 J\y\>i 

It is easy to verify using the definition of G that and assumptions (18) that 
j,,^^eP^G{dx) < CO for all p <E S. and from Theorem 25.17 of Sato [25] this 

means that E[e'''^*] < oo and, moreover, ¥\eP-^^] ~ e**^^^ where 

*("") = oCr^M^+ / (e"^ - I - uxl[^i^i\{x))G{dx) +flu. 

2 ./iB 
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Next, by choosing a > sufficiently large we may ensure that P (supq<j<2- Xt> a) < 
1/2. Then, 

P( inf Xt < -2a) = P( inf Xt < -2a, Xt < -a) 

0<t<T 0<t<T 

+ P( inf Xt < -2a, Xt > -a) 

Q<t<T 

< P(Xt < -a) + P( inf Xt < -2a, Xt > -a), 

Q<t<T 

and for the second term in the preceeding inequality we may use the strong Markov 
property at the stopping time ^ = inf{i > : Xt < —2a} to give P(info<f<TXt < 
—2a, Xt > —a) < l/2P(info<t<T Xt < —2a), and so for any p> Iwe have 

P( inf Xt < -2a) < 2P(Xt < -a) < 2e-P"e'^*(-P). 

0<t<T 



Finally we finish by noting that for e sufficiently small 



( inf St<e) 

0<t<T 



oM.^^^i^n^^"^ 



2eT*(-p) 
qP ' 



Using these facts and the previous theorem we have a random variable tt = iti+tt2, 
where 



Til = 



1 f Wt Jo A^tdt - fo AtdWt J^ Atdt 
aSa 



TloA^dt-[j;^Atdt 



and 7r2 = 7r2.i + 7r2 2 ^I'ith 



7r2,i 



{detCT)^ 



2 I Ajdt I Atdt [ ^r^ 
m \ 2 



Atdt 



dhi 
dhi 



h=0 



Aldt 



h=0 



{A^tfdt 



and 



3^6 



7''2,2 



S-^a 



{detCT) 



Atdt 



9/l2 



T 



Atdt 



9/l2 



/i=0 
T 



A['dt 



h=0 



{A'i)'dt- 



Aidt 







dh2 



h=0 



A't'dt 



such that 



--E[fiST,AT)]^E[f{ST,AT)n]. 



Numerical implementation of these results shows a good degree of accuracy com- 
parable to that achieved by finite difference Monte Carlo in the case of a European 
call. 
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5 Examples 

We show how the formula derived in the previous section should be im- 
plemented to obtain appropriate representations. We will find that the re- 
strictions imposed by Theorems 2 and 3 on the vector fields are often too 
stringent and that we have to get round this problem by localisation. 

5.1 Stochastic volatility models with jumps 

We will consider a volatility process at described by the Heston model 

daf = n{e - crf)dt + rjatdWt. (19) 

We will need the following lemma 

Lemma 1. For any parameter choice with 2k9 > if and for every finite T > 

sup IE[crj" ] < oo. 

0<t<T 

Proof. We let Yt be the squared (5-dimensional Bessel process defined as the 
unique strong solution to the SDE 

t 



Yt = (7l + 6t + 2 I ^/Y.dWs. 
Jo 



For the choice S = ^ we can relate Yt and cr| by the time change (see 
Going-Jaeschke and Yor [19]) 



/ 77^ 



*y ^(e« 

If we let f{t) = f^(e'"* - 1) and T* = f{T) < oo then as t takes values in 
[0, T] so f{t) ranges over [0, T*]. Consequently, 

sup E[a^^] = sup E[(e-'^*y(/(t)))-i] < e^^ sup E[Yt~^]. (20) 

0<t<T 0<t<T 0<t<T* 

Next we notice from the expression for the Laplace transform of Yt (Revuz 
and Yor [23], page 422) 

E[Y-'] = 1^^ E[e-'^^]dX = fi^ + 2Ai)-^/^exp (^^) dX. 

From this and the fact that 6 > 2 we have, for any e > 0, 

sup,<t<T. IE[>"r^] < Ji^il + 2Xt)-^/^dX < (2e)-i. So the proof will be 
complete if we can show 

limsup / (1 + 2Xt)-^/^ cxp -^ dA < oo. 

t-^o Jo \l + 2At/ 
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By using the substitution ^ = 1 — (1 + 2\t) ^ and writing z = t ^ ,y = ^^ 
we need to examine the behaviour of 

^ Jo Jo 

as y ^ oo, and it suffices the check that the expression on the right hand 
side is bounded for large y. To show this, first suppose | — 2 > then we 
trivially have 



Jo 



-2 / (^_^^5/2-2^g-,«^^<^-2(^_^-,)<^- 



Next, suppose | — 2 < 0, then by making the substitution w = (1 — ^) y and 
noticing that for y>\ 



/ (1 - i)^''^''^ye'^^di ^e-y f ', 

Jo Jo 



w'/^-^e^dw 



we see that the right hand side may be bounded uniformly in y since | — 
2>-l. D 

It will be convenient to think of the process at instead, so writing Xt = 
log St to represent the evolution of the logarithm of the stock price the sys- 
tem can be described by the vector SDE 

We shall call this model SVJ. Before the next theorem we introduce the 
notation C^(M'') to indicate the set of real- valued, fc-times differentiable , 
compactly supported functions with domain K'^. We then define I{R'^) to 
be the collection of indicator functions of the form l{a,b): l(a,b]j l[a,b) or l[a,b] 
for some \a\ < |6| < oo and finally a class of real-valued functions on R'', 
JIM''), by 

We will sometimes emphasise the dependence on the initial condition by 
writing Xt = X'^, E = E^, etc.. The following lemma will also be useful 
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Lemma 2. For every y £ R, T > and under the assumption 2k9 > ?; the 
following is true 

lim sup P(X^ e (y - e, y + e)) = 0. 

Proof. We write Jt ~ /g Jg y{^i~v) {dy, ds) and observe that the distribution 
of X^ conditional on Jt and {Wt : < t < T} is Gaussian. Indeed we have 

where 



a= ( r - -a^ ) dt + p atdWt + Jt 

Jo 



This gives 



sup P(X|; e{y-e,y + e))= supE[E[l{x^e(^_,.^+,)j| Jt, {W^t : < t < T}]] 

^ r r^' 1 f {z-a-xf\ ^ 

— sup E / — . cxp dz 

.m [Jy-e V2^/3 V 2/32 ; 

< CeE[/3"i] 

for some constant C < 00. The proof will be complete if we can show 
E[/3^^] < 00, but this is true since the Cauchy-Schwarz inequality gives 



p-' < [T^i^r' 



1/2 



ar^dt 



Then, from the previous lemma. 



E 



1/2- 



ar^dt 



< E 



Jo 



1/2 



< (T sup E[crt-^])i/2 < 00. 

0<t<T 



D 



An application of the extended Bismut-Elworthy-Li formula will give 
the following result. 

Theorem 5. Suppose that the parameters of the SV} model satisfy 2k6 > 77^ and 
/ <E J{M) then, provided \p\ < 1, the following is true 



JrE[/(^T)] = E 



I{St) 



TSo\/l - p^cFs 



-dZ, 
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Remark 7. For the purposes of Monte Carlo applications one would make use of 
the localised Malliavin technique described in Fournie et al [14], and it is clear that 
the class of functions J7'(M) is sufficiently rich for this purpose. In particular, in 
enables us to deal with digital payoffs and European call and put option payoffs. 

Proof Step 1 We assume that / e C^(^) C C'^{^) and note that this implies 
/ e C^ (R) and we let Z? c M be some arbitrary compact subset with x e D. 
It suffices to derive a representation for Xt for / G i7(M), the conclusion for 
St ^vill then follow by applying the result for Xt to the function / o exp e 
J'{M.) , and changing the variable of differentiation to 5*0. 
Step 2 We construct an approximating sequence of SDEs with solution X^ 
such that X^ -^ X a.s. and such that the extended Bismut-Elworthy-Li 
formula can be applied for each X^ . To this end we define for N > 2 



^t''\-{A^ f { --^^(-f) \,.^ rf^i^pH-^) 



^^'"Sr io V3^(-r) - f-f r^ ■ io v o 



ds+ I (V- i^^f v-s ; I j^^ 



'^{PP"f^yW. + fj^[^^]i,-u)[dy.ds) 



(21) 



where the functions h'^ ,g" ,p'^ G C^ (M) are such that 

if \x\ < N 
ii\x\>N + l 




^-^]^ ifx>^ 



and 



p"ix) 



iix< 



4r iix<0 



1 

2N 



where, ^ = ^(f — 1) and, as in Lemma I, S = ^f§-. Moreover, for each A^, 

h^{x) < \x^ for all a; € M, g^ {x) < (^ - ^) ^ for all x £ [0,oo) and 

2^ V x < p^ (x) < 1 for all x € [O, jj] (similar approximating sequences 
for the volatility have been discussed in Ewald [12]). Next, we define the 
stopping times 

TAT = inf I f > : CTt < — I , Cw = inf {t > : ctj > A^} . 

Then, it is well known that for i]^ < 2k9 the volatility never hits zero so 
we have rjv ^ oo a.s. as A^ ^ oo, and since the solution to (19) is non- 
explosive we also have Cw -^ ooa.s. asiV -^ oo. Consequently, for each 
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t G [0, T], X^ = Xt a.s. on the set {t^« > t^N > t} and so X^ -^ Xt a.s. 
as A^ ^ oo. 

Step 3 We confirm that the extended Bismut-Elworthy-Li formula applies 
for each N to deduce 



^E[f{X^)]=E 



r 1 

/(^T ) / ; dZs 



The vector fields driving the SDE defining X^ satisfy the conditions of 
Theorem 2 so we need only verify that the process 

is a martingale for all N. But this is immediate from the fact that the inte- 
grand is bounded (by 2N^/ y^l — p'^). 
Step 4 Next we check that 

To do this we define the sequence of functions (f>^ : D — > M by cf)^ : 
X I— > E^[f{X^)]. We know that each 0^ is differentiable, and it is clear by 
boxinded convergence that (j)^{x) — > (j){x) := E^[/(Xt)] for every x <E D. 
We now confirm that is differentiable with 0'(x) = YnuN ^oo{<j}^ )' {x) ~ 

E[/'(Xt)]. Since, for every N , q^ = 1 this will be achieved if we can 
show 

lim E[/'(X^)]-E[/'(Xt)] (22) 

and the convergence is uniform over x e D. To do this we first show that 
X^ -^ Xt in L^ uniformly in x E D, but since each term E[|Xt — X^W 
is independent of x it suffices the show that X^ -^ Xt in L^, since any 
convergence will then immediately be uniform in x. Before we do this we 
note that a straight forward application of the comparison theorem (page 
269 Rogers and Williams [24]) tells us for each t e [0, T] that yt < a^ a.s. 
where yt is the Ornstein-Uhlenbeck process solving the SDE 

dyt = ~'^ytdt+'^dWt. (23) 

We may also use the proof of the comparison theorem combined with the 
fact that (Tt > a.s. to show that a^ < at a.s. for each t e [0,T]. Then, 
we use the Cauchy-Schwarz and Burkholder-Davis-Gundy inequalities to- 
gether with the fact that [J2i=iXi] — ^ J2i=i ^f ^ show that the family 
{X^ : iV > 2} is bounded in L^. We end up with 

E[{X^y-] < 4x2 + /" ^8j,^2 ^ sTE[h^{aff] + 4E[p^(af )2]) dt 
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Since /i^ (a;) < x^/2,p'^{x) < l + |a;|, supo<t<-r E[(crf)4] < supQ^t^j.E[af] + 
supo<f<y K[yf] < oo we conclude 

E[{X^f]<Ax^ + 8r^T^ + 2T^ sup E[(af )■*] + 8r(l + sup E[{a^f]) 

0<t<T 0<t<T 

and the right hand side of the inequality may be bounded uniformly in N, 
and consequently X^ -^ Xt in L^ uniformly in x. Finally, we verify (22) 
by noting that 



nnx^) .nxT)]\ < Eiwnx^) - /'(Xt)|1{|x«-x.i<.}] 



+ E[\f'{X^)~f{XT)\l^\x-- 



■Xt|>c 



}] 



The first term on the right converging to zero uniformly in x by the uniform 
continuity of /' and the second term likewise by the convergence in proba- 
bility (from Chebyshev's inequality) of X^ to Xt uniformly for x G D and 
the boundedness of /'. 
Step 5 We now establish 



E 



/(^3 



1 



p^(7f 



/o Tv/T 
from which it follows that 
d 



■dZs 



= lim E 



/(^t) 



T^T^p^ia^) 



-dZ,, 



dx 



E[/(Xt)] = E 



f{Xj 



TyT~^ 



-dZ, 



(24) 



To do this we show that K^ —t Kt in L^, where 



K^ = 



1 



, dZs, Kt= , , 

T^/T^pN{a^) Jo T^T^a, 



1 



dZs, 



from which the required result follows immediately by the boundedness 
of /. As a preliminary to this we note that by the time changed squared 
Bessel representation used in Lemma 1 we can deduce 



[tn < T) 



< 



inf a^ < — 

0<t<T N 



inf Yt<—- 

0<t<oo N 



< 



inf r — (e" 
0<t<T \Ak 



< 



N 



iV2€ 



The last line following from the observation that the scale function for a 
(5-dimensional squared Bessel process is s{x) = —x^'^^ (see page 286 of 
Rogers and Williams [24]). We now observe by the Ito-isometry that 



E[\K^~Kt\'] 



1 



< 



T(l-p2 



T(l 

E 



-E 



1 



p^(af) 

2 



P^(^.) 



ds 
1 



p^ias) p^(af). 



ds 
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Using the three facts p'^{x) > x for all x, p'^{(jf) = p^(crt) for t < tn and 
p^{x) > 2^ for all X, we see that 



(l-p2) 

2 



< 



r(i-p2) 

4 



E 



p^(a,) p^(af) 



dn 



{rN<T} 



l-p2 



< 



1 ^ P \0<t<T 



sup ElafT^] + AN^'^FiTN < T) 

0<t<T 

sup EiaiT^] + Ae^^"^] < oo, 



and the fact that K^ -^ Kt in L^ is immediate. 

Step 6 We now relax the regularity conditions on / in two stages. Firstly, we 
extend to / e Cc{^)- To do this we notice that we can identify a sequence 
of functions /„ G C^ (M) with /„ -^ f uniformly and boundedly as n ^ 
oo. The extension is then immediate since bounded convergence implies 
E[/„(Xt)] -^ E[/(Xt)] and, for any compact subset iJ C M, we have 



sup 



d_ 
dx 



nux^)]-n.f{XT)KT] 



< E[if|]i/2 sup E[(/„(X^) - fiX^))']'^' 
xeH 

(25) 



The convergence of the right hand side to zero being immediate from the 
fact that /„ -^ f uniformly, and that X^ = x + S for some random vari- 
able S independent of x. Secondly, we extend to indicator fxmctions of the 
form / = l[a,b\ (the extension to indicators of open and half-open intervals 
being similar). To do this we note that we can construct an approximating 
sequence /„ e Cc{^) having the properties that /„ -^ / pointwise and, 
for any neighbourhoods Ba and Bb of a and h respectively, /„ — / = on 
L := M n -B^ n i?J; for n sufficiently large. We can now repeat the argument 
of the previous paragraph to obtain (25). To show that the right hand side 
of (25) can be made arbitrarily small we let (5 > and fix some e > chosen 
such that 



sup P(X^ e [a - 
xeH 



e, a + e)) < 



4E[X2]i/2 



and 



sup P(X^ e ib - 

x&H 



,b + e))< 



[X|,]l/2 



as we may by Lemma 2. With L = M.r\ {a ~ e, a + e^ {b — e, b + e^ we may 
then choose N such that for all n > A^ we have sup ^^ \fn{y) ~ /(y)| == 
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and we can bound the right hand side of (26) by 



2E[/v^]i/2 g^p (p(^^ e (a - e, a + e)) + F{X^ e{b-e,b + e))) < S. 
xeH 

Since S was arbitrary this completes the result. Since it is clear that (24) 
is stable under taking finite linear combinations the extension to the class 
is immediate. The result for St follows as described in Step 1. D 



Remark 8. By the same argument and under the same conditions as the last 
theroem we can also obtain 



^niisr.. 



-^tJ] = ie 



I{St,. 



,StJ I 
Jo 



Ti 

T^oTT 



1 



9^(78 



-dZ., 



for any n eN andO <Ti <T2 < 



< T„ < T. 



Remark 9. We may apply Theorem 3 together with a similar approximation pro- 
cedure described above to deduce the representation for the gamma 



^niiSr)] = 



E 



(1 - p^)T^Sl 



= E 



1 f^'^ 1 

— dZt j —dZt 

T/2 O"* JO '^t 



ryr 



■ / ^dZAfiSr) 
Jo cr* / 



1 r'^ 1 \ 

^^ 2\rr2Q2 I / — ^^« / — ^^« /('5't) 

(1 - p^)T^S^ \Jt/2 crt Jo <Jt I 



J__d_ 
So OS, 



nfiSr)] 



for f S J'(K). Where, as above, we have initially used Theorem i for Xt and 
deduced the result for St by applying it to the function f o exp and using the 
observation that 






d_ 
'dSo 



d^ 



TT^ — Sq 1^-^ + Sq —^ 



5.2 Stochastic volatility with jumps in the volatility - the 
Matytsin model 

We consider how these ideas may be extended to the model of Matytsin 
where the volatility evolves according to the Heston model with the ex- 
ception that there are jumps which occur in the stock and volatility simul- 
taneously the volatility jumps being of positive deterministic size. This 
volatility process is written as 



dcr* = k( 



- a^)dt + rjatdWt + ^dJt 



where Jt is a Poisson process. With Xt given as in the the SVJ model and a'^ 
as above the pair [Xt^al) describes the Matytsin double jump model ( or 
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SVJJ). Applying Ito's formula we can express the system {Xt, (Jt) in terms 
of our previous notation by the SDEs 



Xt\ f x\ r I ^~¥^s 



<^tj W 'Jo \(f-l)i:--2^s + x{Ja^s-+i-^s-). 



\\ 2 8 J a. 

where £" = M, 7 is the constant jump size in the volatility and /i is a Poisson 
random measure with mean measure v{dy, dt) = XG{dy)dt = Xp{y)dydt 
where here p{y) is the density of the jumps in X. 

Theorem 6. Suppose that the parameters in the SV}} model satisfy 2k0 > tf and 
f S jr(M). Then, provided \p\ < 1, the following is true 



^E[/(5t)] = E 



'^ 1 



f{ST) / -, dZ, 

SoTyJl - p^Os^ 



Proof The proof may be completed by following the steps of the previous 
theorem. The approximating system used is 



Xr\ fx\ . rf r-h^{a^ 



^^ ' - \-oJ " Jo U^(e) - f-f + A(./(af_)2 + 7 - e 



ds 

-^T^- -I- \l . / I rr" 1^ -I- -^z _ /T" II 



;("""<''" ').z.+£(v^''"«)i« 



' f I , — " 

Je 



^j^^r^^_^N_yf^--)idy,ds) 



and then, by the same argument as before, af^ -^ at almost surely. Denot- 
ing a to be the solution of the usual continuous Heston process with the 
same parameters and its approximating process by a^ , and using the fact 
that the jumps in the volatility in Matytsin are non-negative we can apply 
the comparison theorem in between jumps to give the relation 

yt < af < fxf < at (26) 

a.s. for every t e [0,T], where yt is as in (23). Consequently, using {x + 
y)p < 2P^^{xP + yP) we can deduce 

E[ sup (af )^] < 8E[ sup af] + 8E[ sup yf] < 00 
a<t<T o<t<T a<t<T 
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and so a^ -^ at in L^ for every f e [0, T] by dominated convergence. The 
remainder of the argument follows as before, only the elementary observa- 
tion (which follows from (26)) that supQ<(<j.E[cr(~^] < supo<t<7- E[(Tt~^] is 
needed to recycle the estimates already established for the process a to give 
new estimates on the Matytsin volatility a. The extension from C^(]R) to 
J{M) proceeds in the same way as before after Lemma 2 has been verified 
with the Matytsin volatility, which follows from an elementary adaptation 
of the argument given. D 

Remark 10. Under the same assumptions of Remark 8 we may again derive an 
representation for the gamma for the SVJJ model analogous to the one for SVJ. 

6 Numerical Results 

We implement the results for the SVJJ model firstly in the case of a Euro- 
pean call option (payoff {St — K)+) with T = 1 and So = 100 and strike 
K = 100, secondly for a double digital payoff of the form liKi.K2]iST) 
again with T ~ 1,5*0 = 100 and Ki ~ 100, K2 = HO. Finally, we implement 
for the delta of a digital Cliquet option with payoff profile 1[k'.k:;]{St — 
Sti), where T = l,Ti = 0.5, K^ = 5 and A'^ = 10. The model parameters 
we use are r = 0, p = -0.7, 7 = 0.4, erg = 0.1, A = 1.0, 9 = 0.08, k = 4.0, 
77 = 0.6, and we assume that the jumps in logS* are distributed normally 
with mean —0.1 and standard deviation 0.1. 
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Figure 1: Gamma for a European call option with parameters as above 
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Figure 2: Delta for a double digital option with parameters as above 
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Figure 3: Delta for a digital Cliquet option with parameters as above 
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